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Together, the standard theory of big-bang nucleosynthesis (BBN) and the primeval deu- 
terium abundance now very precisely peg the baryon density. Based upon our analysis of 
the deuterium data and the theoretical uncertainties associated with the BBN predictions, 
we determine CIbH 2 — 0.020 ± 0.002 (95% C.L.), with the uncertainty from the measured 
deuterium abundance about twice that from the predicted abundance. We discuss critically 
the reliability of the BBN baryon density, and in light of possible systematic uncertainties 
also derive a very conservative range. We conclude that within the standard cosmology and 
standard theory of BBN a baryon density Q,sh 2 — 0.032 (the central value implied by recent 
CMB anisotropy measurements) simply cannot be accommodated. 
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I. INTRODUCTION 

Since the determination of its solar system abun- 
dance by the Apollo astronauts almost thirty years 
ago, deuterium has been used to constrain the den- 
sity of ordinary matter [QJ^]. The reason is simple: 
big-bang deuterium production has a strong depen- 
dence upon the baryon density (oc Pg 1 ' 6 ), and as- 
trophysical processes since have only destroyed deu- 
terium 

Until recently, deuterium was used to set an upper 
limit to the baryon density (around 10% of critical 
density), based upon the fact that the big bang must 
produce at least the amount of deuterium seen in 
the local interstellar medium (ISM). Together with 
measurements of the abundances of the other light 
elements produced in the big bang, a concordance 
interval for the baryon density was derived, 0.007 < 
il B h 2 < 0.024 B. For two decades the big-bang 
nucleosynthesis (BBN) baryon density has stood as 
the best determination of the amount of ordinary 
matter and the linchpin in the case for nonbaryonic 
dark matter || . 



A dramatic change occurred in 1998 when the 
abundance of deuterium was measured in high- 
redshift clouds of pristine gas backlit by even more 
distant quasars. Tytler and his collaborators have 
now determined deuterium abundances for three 
clouds and derived upper limits for a number of 
other clouds Based upon these and other 

results we infer a primeval deuterium abundance, 
(D/H) P = (3.0±0.4) x 10~ 5 (95% C.L.), which leads 
to a precision determination of the baryon density 

1: 

n B h 2 = 0.020 ±0.002 (95% C.L.). (1) 

£Ib is the fraction of critical density contributed 
by baryons, h = i/o/100kmsec~ 1 Mpc~ 1 , and 
the physical baryon density ps = 1.878(r2s/i 2 ) x 
10- 29 gcur 3 . 

A key test of the BBN prediction, and indeed 
the consistency of the standard cosmology itself, lies 
ahead. Measurements of cosmic microwave back- 
ground (CMB) anisotropy on small angular scales 
can ultimately determine the baryon density to an 
accuracy of around 1%. The physics involved is 
very different: gravity-driven acoustic oscillations 



of the photon-baryon fluid when the Universe was 
around 500000 years old. The first step toward 
this important goal was taken recently when the 
BOOMERANG and MAXIMA CMB experiments 
reported results for the baryon density: f^/i 2 = 



0.032 



+0.009 
-0.008 



(95% C.L.) 



At about the 2a level, 



this independent measure of the amount of ordinary 
matter agrees with the BBN prediction and supports 
the longstanding BBN argument for non-baryonic 
dark matter. 

The CMB determination of the baryon density 
should improve dramatically over the next few years, 
making a very precise comparison of the two meth- 
ods possible. The difference between the BBN and 
CMB baryon densities has already triggered lively 
discussion in the literature p}]. Motivated by this, 
we have written this paper to explain the BBN pre- 
diction for the baryon density and the associated 
uncertainties. 
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FIG. 1. Cross-section data and fits for one of the pro- 
cesses important for determining the BBN deuterium 
yield. Solid curves indicate best-fit and 95% C.L. lim- 
its from our Monte Carlo method |jll|jl^| ; broken curves 
indicate the corres pon ding best-fit and 95% C.L. curves 
from earlier work [n3|. The bracket at the bottom indi- 
cates the energy range where this cross section is needed 
in order to compute all light-element abundances to an 
accuracy of one-tenth of their current uncertainties. 



II. UNCERTAINTIES IN THE PREDICTED 
DEUTERIUM ABUNDANCE 

The standard scenario for big-bang nucleosyn- 
thesis begins with the assumptions of (i) the 
isotropic and homogeneous Friedmann-Robertson- 
Walker (FRW) cosmology, (ii) three massless (or 
very light) neutrino species, (iii) zero (or very small) 
chemical potentials for the neutrino species, and (iv) 
spatial homogeneity of the baryon density. Within 
the standard theory, the predictions for the light- 
element abundances depend only upon the baryon 
density (more precisely, baryon-to-photon ratio) and 
the dozen nuclear-reaction rates that enter the cal- 
culation. 

Essentially all of these rates have been measured 
at energies relevant for BBN, and because large num- 
bers of measurements exist for most of the cross sec- 
tions, reliable estimates of the uncertainties can be 
made (see Fig. [I]). We have used all the extant nu- 
clear data to estimate directly the uncertainties in 
the BBN predictions 0,0]. We integrated the nu- 
clear data over thermal distributions to obtain ther- 
mally averaged rates. Then, uncertainties in the pre- 
dictions were estimated by a Monte Carlo method, 
repeatedly calculating the abundances with individ- 
ual cross-section data being varied according to their 
quoted experimental uncertainties (also accounting 
for correlated normalization errors). From the re- 
sulting distribution of predicted abundances we have 
derived the uncertainties. 
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FIG. 2. D/H uncertainties (from left to right): mea- 
sured D/H 0, uncertainty estimated by SKM Jl3| , 
our Monte Carlo uncertainty, and uncertainties arising 
from the individual reactions, all at 95% C.L. and for 
Q B h 2 = 0.020. 
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The uncertainty in the D/H prediction around 
Fish 2 = 0.020 arises mainly from the reactions 
d(d, n) 3 He, d(p, 7) 3 He, and p(n,j)d, with a smaller 
contribution from d(d,p)t (see Fig. ||). The total un- 
certainty from all of these combined is 5.4% at 95% 
C.L., and increases slightly for larger Og/i 2 . Com- 
bining this with the slope of the power-law depen- 
dence of D /H on baryon density, we obtain a nuclear 
contribution to the uncertainty in the deuterium- 
derived baryon density of 3.4% (95% C.L.). 

We end this discussion with two comments. First, 
BBN actually determines the baryon-to-photon ra- 
tio when the Universe was about 100 seconds old, 
r] = (5.6 ±0.6) x 10~ 10 (95% C.L.); to translate this 
into a baryon density today two things are needed: 
(i) average mass per baryon; and (ii) the assump- 
tion of adiabatic expansion since BBN (l4j . At the 
1% level, the average mass per baryon depends upon 
chemical composition; for the post-BBN primordial 
mix or universal solar abundance, m = 1.670 — 
1.671 x 10~ 24 g. Adiabatic expansion since BBN is 
a feature of the standard cosmology. With these 
assumptions, fl B h 2 = (3.650 ± 0.008) x 10 7 r? (95% 
C.L.), where the error comes from the uncertainty in 
the CMB temperature, T = (2.725 ± 0.002) K (95% 
C.L.) @. 

The second comment is an explanation of why 
our estimates of the light-element uncertainties are 
about a factor of two to three smaller than the pre- 
vious very thorough study by Smith et al. (SKM) 
p3| in 1993. Some of it is simply improved measure- 
ments. However, the bulk of the difference involves 
technique. As described, our analysis used all the 
data directly, each data point weighted by its error 
bar. SKM estimated 2ct limits to cross sections by 
constructing envelopes (motivated by theory) forced 
to contain most of the measurements, and then from 
these, derived overall rate uncertainties. Although 
much simpler to implement, this technique gives too 
much weight to experiments with large error bars 
and to cross-section data in energy intervals that 
are not as important to BBN, and results in unnec- 
essarily conservative error estimates. 

III. UNCERTAINTIES IN THE OBSERVED 
DEUTERIUM ABUNDANCE 

In 1976 Adams jl6| pointed out that the primor- 
dial abundance of deuterium could be determined by 
observations of Lyman-series absorption of quasar 
light by intervening high-redshift (z > 2.5) clouds 
of pristine gas, by using the fact that the deuterium 



feature is isotopically shifted by 82km/s to the blue 
(see Fig. ||). Successfully implementing this idea 
awaited the advent of the 10-meter Keck telescope 
and its HIRES spectrograph jl7|] . At present, four 
absorbing gas clouds place stringent constraints on 
the primeval deuterium abundance, and another five 
provide independent consistency checks. Here we 
briefly summarize the state of the observations, and 
discuss possible systematics. For a detailed review, 
see Ref. M. 

The technique is simple and straightforward; the 
relative abundances of deuterium and hydrogen are 
measured through absorption profile fitting of the 
Lyman series. In these clouds of virtually primordial 
material (metal abundance less than 1% of solar), 
the relative column densities of neutral deuterium 
and neutral hydrogen yield the primeval D/H ratio 
without correction for ionization or destruction of 
deuterium by stars ^7j. The major obstacle is the 
discovery of systems that are suitable for deuterium 
detection. To observe the weak deuterium feature, 
the neutral hydrogen column density must exceed 
about 10 17 cm 2 , and such clouds are relatively rare. 
The current rate of success is approximately 4 suit- 
able clouds per 100 quasars studied. (Note, for each 
high-redshift quasar there are hundreds of interven- 
ing gas clouds, but typically only one with suffi- 
ciently high column density to see deuterium.) 

There are three systems at z > 2.5 where deu- 
terium has been detected: D/H = (3.3 ±0.6) x 10~ 5 
PJ; (4.0 ± 1.4) x 10- 5 @ and 2.5 ± 0.5 § (95% 
C.L.). The third and newest system has a much 
higher hydrogen column density, and deuterium is 
seen in Lyman-/?, 7, 5, 6, and 7. A fourth system 
provides a strong upper limit: D/H < 6.7 x 10~ 5 
(95% C.L.) gfj. 

Using the likelihood distributions for the first two 
detections and a gaussian likelihood for the new de- 
tection, we infer (D/H) P = (3.0 ± 0.4) x 10~ 5 (95% 
C.L.). We note that the dispersion of the three de- 
tections is somewhat larger than one might expect 
based upon the estimated error of the mean, though 
not unreasonably so: The reduced x 2 = 7.1 for two 
degrees of freedom has a probability of 3%. 

Another five extragalactic absorbers that give lim- 
its on D/H at high redshift have been analyzed 
[^l],^2). They add no weight to the D/H determina- 
tion, but are all consistent with (D/H)p = 3 x 10~ 5 . 
Fig. ^ summarizes the deuterium detections and up- 
per limits. 
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FIG. 3. Spectrum of Q1937-1009; blueward of the 
characteristic Lyman-a emission line of the quasar is the 
"forest" of Lyman-a absorption due to the hundreds of 
intervening gas clouds. The lower panel shows a blowup 
of the region around the deuterium detection, a cloud at 
redshift z = 3.572, and the model fit. 

Earlier claims of a factor of ten higher deuterium 
abundance in an absorption cloud associated with 
Q0014+8118 (23| were shown to be a result of 
misidentification of the putative deuterium feature 
p4j . The aforementioned upper-limit systems also 
argue strongly against high deuterium. In addition, 
recent analyses of a low redshift system suggestive 
of high D/H [^5 26 are unconvincing due to lack of 
spectral coverage of the entire Lyman series. 

Although the technique is simple and direct, there 
are two important sources of systematic uncertainty. 
The first arises from "hydrogen interlopers," low 
column-density gas clouds that are coincidentally 
situated to mimic deuterium absorption |^2],^6] p8| . 
Interlopers result in a one-sided systematic bias 
which can lead to an overestimation of D/H. Ab- 
sorbing clouds are ubiquitous at high redshifts 
and the chance probability of an interloper is non- 
negligible. 

The probability of an interloper depends on many 
factors, including redshift, neutral hydrogen column 



density, and velocity dispersion as well as the intrin- 
sic value of D/H. Here, we present an a posteriori 
estimate of the probability for hydrogen contamina- 
tion at the level of 10% or more in each system. In 
the three absorption systems with detections, the a 
priori probability for 10% contamination is: 0.007, 
0.01 and 0.001 respectively, assuming that the inter- 
loper must fall within 5km/s of the expected posi- 
tion of deuterium absorption. (Because of its very 
high hydrogen column density, the third system is 
particularly immune to the possibility of an inter- 
loper.) Multiplying these a priori probabilities by 
the number of clouds searched (about 25), we ar- 
rive at estimates for the a posteriori probabilities, 
0.18, 0.25, 0.03 respectively. The final joint prob- 
ability that all three systems are contaminated by 
more than 10% is less than 0.1%. 

This small probability, and the consistency of the 
three deuterium detections, argues strongly for deu- 
terium detections rather than hydrogen interlopers. 
One can attempt to correct for hydrogen contamina- 
tion by assuming a flat prior for contamination of 
to 10% in the lower column-density systems, Q1009 
and Q1937. The likelihood for (D/H)p, in this case, 
has a central value that is 3% lower and an uncer- 
tainty that is slightly (20%) larger. This analysis 
further suggests that contamination by interlopers 
is not a significant systematic; further, allowing for 
this small contamination lessens the possible signif- 
icance of the larger than expected scatter in the de- 
rived deuterium abundances. 
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FIG. 4. Summary of deuterium detections and upper 
limits as a function of the hydrogen column density (all 
at 95% C.L.). The diagonal line indicates the minimum 
deuterium abundance that could be detected as a func- 
tion of hydrogen column density. Note that the newest 
system lies well above the detection threshold. 
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The other important source of systematic error 
involves model fitting of absorption line profiles in 
quasar spectra. A small number of Voigt profiles 
were used to model the deuterium absorption sys- 
tems jl^,[l{|. The effect of model-parameter choice 
has been studied extensively in the measurement of 
D /H at high redshi ft |22| 30| ] , as well as in the local in- 
terstellar medium |31fl3S "including the uncertainty 
of the unabsorbed quasar continuum, the estimated 
possible systematic bias is at the level of 5%. This 
bias could affect all sightlines, but it should have a 
different magnitude and direction from one system 
to the next. 

Finally, it has been argued that astrophysical 
mechanisms that significantly alter the deuterium 
abundance in these clouds are implausible p7| ; the 
relatively small intrinsic scatter allowed by the cur- 
rent data argues against any significant post-BBN 
production or destruction, as more scatter would 
be expected. This argument is weaker than those 
above, since there could be some real scatter in the 
data. 



IV. WHY BBN CANNOT ACCOMMODATE 

Q B h 2 = 0.03 



processes since BBN have only destroyed deuterium. 
This is in stark contrast to the abundance of heavy 
elements (around 2%) which indicates about half the 
material has been processed through stars, and thus 
implies a primeval deuterium abundance of about 
3 x 1CT 5 . Further, the D/H prediction at n B h 2 = 
0.03 is also lower than the abundance inferred for 
the pre-solar nebula, (D/H)= (2.1±1.0) x 10~ 5 (95% 
C.L.) @. 

The predicted 7 Li abundance is about six times 
that measured in the atmospheres of old pop II halo 
stars, ( 7 Li/H)= [1.7 ± 0.1 (95%C.L.) ± 0.2 (syst)] x 
10~ 10 [f||. While there is lively debate about the 
possible depletion of Li in these stars (by rotation- 
ally driven mixing and convective burning), there is 
a consensus that any such dep letion must be less 
than about a factor of two |3(|. Thus, 7 Li too is a 
serious problem for this baryon density. 

Finally, though there is no consensus about the 
primeval abundance of 4 He, and discussion contin- 
ues about possible systematic error, Yp = 0.250 is 
inconsistent with the two largest compilations of He 
measurements. The primordial 4 He abundances in- 
ferred from these two studies of metal-poor, extra- 
galactic Hn regions are: Y P = 0.244 ±0.004 Q and 
Y P = 0.234 ± 0.004 M, both at 95% C.L. 



Because the uncertainty associated with the CMB 
baryon density is significant (about 15%), there is 
nothing special about the central value, tt B h 2 = 
0.032. Nonetheless, it is interesting and instruc- 
tive to ask if BBN could accommodate such a high 
value. The answer is a resounding no for three 
reasons: D, 7 Li and 4 He, in that order. For a 
baryon density of VL B h 2 = 0.03, the predicted abun- 
dances are: (D/H) = (1.6 ± 0.12) x lO" 5 ; ( 7 Li/H) 
= (8.5 ± 0.7) x 10- 10 ; and Y P = 0.251 ± 0.001 (95% 
C.L.). All three conflict significantly with observa- 
tions. 

Aside from high-redshift hydrogen clouds, the 
most precise determination of the deuterium abun- 
dance is that measured with the HST for the local 
ISM, (D/H)= (1.5±0.2) x 10~ 5 (95% C.L.) Q. The 
technique is the same, except on a smaller scale: ab- 
sorption by clouds of neutral gas along the lines of 
sight to many nearby stars is measured. (Variations 
in the local D/H, both upward and downward, have 
been reported and may be statistically significant 
p2| , p3]| . If real, they probably reflect the inhomo- 
geneity of the local ISM.) The ISM value for D/H is 
essentially equal to the predicted primordial abun- 
dance for Q, B h 2 = 0.03 and would imply that the 
local ISM is pristine material, because astrophysical 



V. A VERY CONSERVATIVE RANGE FOR 
THE BBN BARYON DENSITY 

We believe our stated range for the BBN baryon 
density, tt B h 2 = 0.020 ± 0.002 (95% C.L.) is well 
justified. However, two recent reviews have quoted 
broader 95% C.L. intervals: 0.015 < n B h 2 < 0.023 



59| and 0.004 < il B h 2 < 0.021 |40|. The differ 



ences are simple to explain: The first review |39[ 
used the older SKM analysis for the error in the 
predicted deuterium abundance (see Fig. ||). The 
second review |^(| still allows for the possibility that 
the primeval D/H is as large as 3 x 10 -4 . As de- 
scribed above, we believe the case for high primeval 
D/H is simply no longer tenable. 

The question remains, just how high can the 
baryon density be pushed. To answer that ques- 
tion and obtain a "very conservative range" for the 
baryon density, we revert to the older, more con- 
servative SKM analysis of the input nuclear data - 
though unlikely, our Monte Carlo analysis could be 
driven by an experiment with understated errors - 
and use a very conservative range for the primeval 
deuterium abundance, (D/H) P = (2.1 - 3.9) x 10~ 5 . 

This range for the primeval deuterium abundance 
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is derived from a weighted average of the three de- 
tections, with the error of the mean being estimated 
from the standard deviation of the three detections 
using the standard formula. By so doing, we are 
in effect disregarding the estimated errors for the 
individual measurements, and instead using the dis- 
persion of the measured values to estimate the error 
of the mean. If the statistical errors are larger than 
estimated or if there are unknown systematic errors, 
this approach might more accurately reflect the un- 
derlying errors. (We note that O'Meara et al. [0, 
motivated by the larger than expected dispersion of 
the three detections, have advocated a similar ap- 
proach for estimating the uncertainty in the primeval 
deuterium abundance.) 

Adding the theoretical and observational uncer- 
tainties discussed above in quadrature, we obtain a 
very conservative range for the BBN baryon den- 
sity: fish 2 = 0.016 — 0.026, which is about twice 
the width of what we believe to be a well justi- 
fied 95% confidence interval. Precisely because this 
very conservative range is predicated upon the pos- 
sibility of systematic error, or error that is not well 
quantified by a Gaussian distribution function, we 
have not assigned a confidence level to it. We also 
note that at the upper limit of this extreme range, 
rish 2 — 0.026, all the light-element abundances, 
(D/H) P = 2.0 x 10~ 5 , ( 7 Li/H)= 6.6 x lO" 10 and 
Yp = 0.249, are uncomfortably different from their 
inferred values. 



VI. CONCLUDING REMARKS 

A flood of high-precision data - from CMB mea- 
surements to very large redshift surveys to preci- 
sion D/H determinations - is transforming cosmol- 
ogy. Soon, no one will remember a time when the 
phrase precision cosmology was an oxymoron. The 
new data are testing our most promising ideas about 
the early Universe as well as the consistency of the 
big-bang framework itself. A comparison between 
the BBN and CMB baryon densities will be one of 
the most important consistency tests. 

The use of BBN to determine the baryon density 
is a mature subject. The predictions have been scru- 
tinized and the required nuclear data are measured 
at the relevant energies. The determination of the 
deuterium abundance in nearly pristine high-redshift 
gas clouds now involves nine systems, and consid- 
erable attention has been paid to systematic error. 
While our stated uncertainty in the BBN baryon 
density is small, fl B h 2 = 0.020±0.002 (95% C.L.), it 



is well justified (our very conservative range is only 
twice as broad) . As more deuterium systems are dis- 
covered and analyzed, the uncertainty will shrink. 
Even without improvement in the nuclear data, a 
precision of about 4% (at 95% C.L.) is possible. 

On the other hand, CMB anisotropy measure- 
ments have just recently achieved sufficient angular 
resolution to probe the baryon density. The first 
result, a 15% determination of the baryon density, 
fish 2 = 0.032t^oo8' i s encouraging. It supports the 
longstanding BBN argument for nonbaryonic dark 
matter and agrees with the BBN baryon density at 
about the 2a level. It is important to note that 
the CMB results for Uph 2 depend upon the number 
of free cosmological parameters used in the analysis 
and the priors assumed for them 0. The CMB de- 
termination will improve significantly over the next 
few years as more experiments probe the sub-degree 
angular scales, making a more robust comparison 
with the BBN prediction possible. 

As we have emphasized, there are three good rea- 
sons why BBN cannot tolerate VLsh 2 — 0.03: the 
predicted abundances of D, 7 Li and 4 He; all would 
conflict significantly with observed abundances. The 
standard cosmology and standard BBN cannot - and 
probably will not have to - accommodate a baryon 
density this large. 

Should future CMB measurements zero in on a 
baryon density higher than our extreme upper limit, 
fish 2 = 0.026, they would, in our opinion, implicate 
nonstandard cosmology or BBN. Among the possi- 
bilities are: large neutrino chemical potentials [ ElJ , 
entropy reduction since BBN due to exotic physics 
Jl4[ , a decaying tau neutrino [Q, neutrino oscilla- 
tions (i^] , or the inconsistency of the standard cos- 
mology. At the moment, none of these possibilities 
are particularly compelling, and we will wait to see 
if they are necessary. 
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